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Abstract 



In this chapter first the statistics of the standard and truncated Pareto distributions are de- 
rived and used to fit empirical values of asteroids diameters from different families, namely, 
Koronis , Eos and Themis, and from the Astorb database. A theoretical analysis is then car- 
ried out and two possible physical mechanisms are suggested that account for Pareto tails 
in distributions of asteroids diameter. 
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1. Introduction 

The interest in the study of asteroids in the inner solar system lays in their connection with 
the formation of planets and their temporal evolution. 

Among others, studies on asteroids formations and evolutions involve 

1. the effects of asteroid collisional history on sizes and spins of present-day objects HI, 

2. realistic collisional scaling laws and the implications of including observables, such 
as collisional produced families, in constraining the collisional history of main-belt 
asteroids (H, 

3. the creation of a model of the main asteroid belt whose purpose is to describe distri- 
bution of size size frequency of asteroids and simulate their number pQ. 

4. temporal evolution for 4.2 Myr of test particles, which were initially placed on a 
perfectly rectangular grid and subjected to gravitational interactions with the Sun and 
five planets, from Mars to Neptune , see EJ. 

On the other hand, it has been shown that experimental observations can be fitted with 
a differential size distribution 

n = dN/dD = n(D) oc D~ a , (1) 

where D is the diameter in Km , a the exponent of the inverse power law and n the num- 
ber of asteroids comprised between D and D + dD. Measurements of the properties of 
13,000 asteroids detected by the Sloan Digital Sky Survey (SDSS) present a differential 
size distribution that for D > 5Km is n oc I)" 4 and for D < 5Km is n oc D , see [4J. 

The ongoing simulations as well the observations require a careful analysis of the Pareto 
distribution 0|6l and the truncated Pareto distribution 0[U[9]]. This paper presents in Sec- 
tion^ a comparison between the Pareto and the truncated Pareto distributions. In Sectionl3~1 
the theoretical results are applied to the distribution of the radius of asteroids. Two phys- 
ical mechanisms that produces a Pareto type distribution for diameters are presented in 
Section 0] 



2. Statistical Distribution 

Let X be a random variable taking values x in the interval [a, oo], a > 0. The probability 
density function (in the following pdf) named Pareto is defined by f6\ 

f{x;a,c) = ca c X - {c+l) , (2) 

c > 0, and the Pareto distribution functions is F(x : a, c) = 1 — a c x~ c 

An upper truncated Pareto random variable is defined in the interval [a, b] and the cor- 
responding pdf is 

coPx ( c ~^^ 

f T (x;a,b,c) = _ (a . c , (3) 
1 Kb) 
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[9] and the truncated Pareto distribution function is 

F T (x;a,b,c) = — ^ . (4) 

Momenta of the truncated distributions exist for all c > 0. For instance, the mean of 
fr(x; a, b, c) is, for c / 1 and c = 1, respectively, 

ca l-il) c - 1 ca c , b 



Similarly, if c ^ 2, the variance is given by 

ca 2 1- 
(c-2) 1 



2 ca 2 l-(f) c - 2 2 
* = 7 ^ , L < S > 2 , (6) 



whereas for c = 2 

T - T|f ln--<x> 2 . (7) 
In general the n — th central moment is 

r-b 



(x- < x >) n f T {x)dx 



(- < x >) n a- c 2 F 1 (-c, -n; 1 - c; -^-^) ((a 6 )" 1 - (ft 6 )" 1 ) 
-(-<x>)"6- c 2 F 1 (- C ,-n; 1 - c; — ) ((a )" 1 - (ft")" 1 ) ~* (8) 



where 2-^1(0, b; c; 2) is a regularized hypergeometric function, see [fT0lfTTHT2l . An analo- 
gous formula based on some of the properties of the incomplete beta function, see lfT3l and 
fffl , can be found in 03). The median m of the Pareto distribution is 

m = 2 c ' 1 a , (9) 

and the median of truncated Pareto mr 

m T = a2 c ~ 1 (a c b- c + l)~ c_1 . (10) 



Parameters of the truncated Pareto pdf from empirical data can be obtained via the 
maximum likelihood method; explicit formulas for maximum likelihood estimators (MLE) 
are given in Q, and for the more general case in 0, whose results we report here for 
completeness. 

Consider a random sample X = x\, x 2 , ■ ■ ■ , x n and let xn\ > xr 2 ) > • • • > x (n) denote 
their order statistics so that xm = max(xi, x 2 , • • • , x n ), xt n ) = min(a?i, x 2 , ■ ■ ■ , x n ). 
The MLE of the parameters a and b are 

a = X(n), b = x {1) , (11) 
respectively, and c is the solution of the equation 
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i + rz ?m # - 5> x * - lnx w ] = °' (12) 
'nil' * =i 

& 

There exists a simple test to see whether a Pareto model is appropriate [9]: the null 
hypothesis Hq : v = oo is rejected if and only if xh\ < [nC/(— In q)] 1 ^, < q < 1, 

where C = a c . The approximate p-value of this test is given by p = exp j— nCx^j, and 
a small value of p indicates that the Pareto model is not a good fit. 

Given a set of data is often difficult to decide if they agree more closely with / or 
fx, since, in the interval [a, b], they differ only or a multiplicative factor 1 — (a/b) c , that 
approaches 1 even for relatively small values of c if the interval [a, b] is not too small. For 
this reason, rather than / and fx, the distributions P(X > x) and Pr(X > x) are used, 
often called survival functions, that are given respectively by 

P(X >x) = S(x) = 1 - F(x; a, c) = a c 2T c (13) 

and 

P T (X >x) = S T (x) = 1 - F T (x; a, b, c) = — {x ~] \~ .f^ • (14) 

1 \b) 

The Pareto variate X can be generated by 

X : a,c~ a(l - R)~° , (15) 
and the truncated Pareto variate Xt by 

X T :aAc~a(l-R(l-(^) c )^ C , (16) 
where R is the unit rectangular variate. 



3. Application to the asteroids 

We have tested the hypothesis that diameters of asteroids follows a Pareto distribution by 
by considering different families of asteroids, namely, Koronis , Eos and Themis. 

The sample parameter of the families are reported in Table [TjTable [2] and Table [3] , 
whereas Figure [H Figure [2j Figure [3] report the graphical display of data and the fitting 
distributions. 

Table 1. Coefficients of diameter distribution of the Koronis family . The parameter c is 
derived through MLE and p=0.033 . 



a [km] 


b [km] 


c 


n 


P(X > x) 


25.1 


44.3 


3.77 


29 


truncated Pareto 


25.1 


oo 


5.04 


29 


Pareto 
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In (D(km)) 

Figure 1. ln-ln plot of the survival function of diameter distribution of the Koronis Family: 
data (empty circles), survival function of the truncated Pareto pdf (full line) and survival 
function of the Pareto pdf (dotted line). A complete sample is considered with parameters 
as in Table [TJ 




In (D(km)) 

Figure 2. ln-ln plot of the survival function of diameter distribution of the Eos Family: 
data (empty circles), survival function of the truncated Pareto pdf (full line) and survival 
function of the Pareto pdf (dotted line). A complete sample is considered with parameters 
as in Table 12 
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Table 2. Coefficients of diameter distribution of the Eos family . The parameter c is derived 
through MLE and p=0.681 . 
a [km] b [km] c n P(X > x) 

30.1 flO 3^80 53 truncated Pareto 

30.1 oo 3M 53 Pareto 



Table 3. Coefficients of diameter distribution of the Themis family . The parameter c is 
derived through MLE and p-0.67 . 
a [km] b [km] c n P(X > x) 

35.3 249 Z5 53 truncated Pareto 

35.3 oo 2.6 53 Pareto 



In case of the Koronis family Pt fits the data better than P and indeed p = 0.039 
is correspondingly small, whereas for the Eos family, P performs slightly better than Pt 
(p=0.68), and the estimated of c are very closed in both cases. Finally in the third case, the 
Themis family, the two distributions are the same, due to the fact that the ratio a/6 = 0.14 
is small. 

Another interesting catalog is the Asteroid Orbital Elements Database (Astorb) which 
is visible at http://vizier.u-strasbg.fr/; the sample parameter of the asteroids with diameter 
> 90 Km is reported in Tableland the fitting survival function in Figure 0] 

Table 4. Coefficients of diameter distribution of the asteroids extracted from Astorb 
database with diameter > 90 Km. The parameter c is derived through MLE and p=0.53 . 



a [km] 


b [km] 


c 


n 


P(X > x) 


90.59 


848.4 


2.71 


272 


truncated Pareto 


90.59 


oo 


2.75 


272 


Pareto 



4. Simulating Pareto tails 

Two models that explains the Pareto tails are presented. The first analyzes the possibility 
that the asteroids are formed from smaller bodies and the second one introduces a fragmen- 
tation model at the light of the Voronoi diagrams. 

4.1. Accretion 

As a simple example of how a distribution with power can be generated, consider the growth 
of a primeval nebula via accretion, that is the process by which nebulae "capture" mass. We 
start by considering an uniform pdf for the initial mass of N primeval nebulae, m, in a range 
m-min < m < m max . At each interaction the i-th nebula has a probability Aj to increase 
its mass m^ that is given by 



Aj = (1 — exp(— akrrii)), 



(17) 
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Figure 3. ln-ln plot of the survival function of diameter distribution of the Themis Family: 
data (empty circles), survival function of the truncated Pareto pdf (full line) and survival 
function of the Pareto pdf (dotted line). A complete sample is considered with parameters 
as in Table [3] 




n (D(km)) 



Figure 4. ln-ln plot of the survival function of diameter distribution of the asteroids 
extracted from Astorb database with diameter > 90 Km, data (empty circles), survival 
function of the truncated Pareto pdf (full line) and survival function of the Pareto pdf (dotted 
line). A complete sample is considered with parameters as in Table |4] 
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In (D(km)) 



Figure 5. log-log plot of the survival function of the diameter distribution for the primeval 
nebula. The truncated Pareto parameters are c=2.75 and p=0.00026 . 

where ak is a parameter of the simulation; thus more "massive" nebulae are more likely to 
grow, via accretion. The quantity of which the primeval nebula can grow varies with time, 
to take into account that the total mass available is limited, 



where <5m(0) represents the maximum mass of exchange and r the scaling time of the 
phenomena. The simulation proceeds as follows: a number r, is randomly chosen in the 
interval [0, 1] for each nebula, and, if r < Aj, the mass m, is increased by 5m(t), where 
t denotes the iteration of the process. The process proceed in parallel : at each temporal 
iteration all the primeval nebulae are considered. 

Results of the simulations have been fitted with both Pareto survival distributions, see 
Figure [5] 

4.2. Fragmentation 

The distribution of fragments size arising from breaking of material is still subject of re- 
search as it depends on the actual fragmentation process. The first model was developed 
in lfl6l : it assumes a time dependent probability of fracture of a ring after a critical strain 
has been reached in the material. The resulting distribution show that the frequency of 
occurrence of fragments masses follows a cumulative distribution of the form : 



5m(t) = 5m(0) exp(— tjr) 



(18) 




(19) 



N 



where N m is the number of fragments each of whose mass is greater or equal to m , N 
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is the total number of fragments , 



m 

« = y , (20) 



and m is the averaged mass of a fragment. 
4.2.1. Fractal distribution of fragments size 

In order to generate cells resulting in a fractal distribution of their volumes the following 
method can be adopted, which generalizes the procedure presented in [17]. 

Consider a domain V subdivided into iVo cubic cells with a linear dimension I that, that 
in the following, will be referred as zero-order cells. A zero-order cell is divided into A; 3 
smaller cubes called zero-order elements, with linear dimension / / k and volumes 

V Q 

Vi = (21) 

where Vq is the volume of zero-order cell. UP is the probability of a zero-order cell to 
be fragmented, the number N\ of zero-order elements generated by fragmentation is given 
by 

N\ = P • k 3 N , (22) 
and the number No a of zero-order cells that have been not fragmented is 

N 0a = (1 - P)N . (23) 

Each zero-order element now becomes a first-order cell that can be fragmented into 
first-order elements of volumes 

« = w? <24) 

and the number of fragmented first-order elements is 

N 2 = Pk^Nx = (k 3 P) 2 N . (25) 
The number iVi a of first order cells that have not been fragmented is given by 

iV la = k 3 P(l - P)N . (26) 

Now first-order elements can be considered second-order cells and the procedure re- 
peats itself. The volume of the nth-order cell V n is 

and, after fragmentation, the number of nth-order cells N na is 

N n = (Pk 3 ) n N 0a = (Pk 3 ) n (l - P)N (28) 
Taking the natural logarithm Eqs. (T27T ) and (T28T ) leads to 

In— = -nln(A; 3 ), (29) 
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In (D(km)) 

Figure 6. log-log plot of the survival function of the diameter distribution as given by the 
Voronoi Diagrams in presence of 1000 fractal seeds. The parameters of the simulation are 
k =2 , P = 0.92 and therefore Dj r =2.80. The truncated Pareto parameters are c=1.52 and 
p = 0. 



ln na = _ nln (p fc )3 _ (30) 

^From Eqs. (|30l and ( f30b it is straightforward to obtain, by elimination of n : 

ln[Pfc 3 ] 



N nn . X, 



W*?] , (31) 



N 0a Vq 

that is a fractal distribution with dimension D given by 

D =n^Wr ■ <32) 

Now we can consider the center of cells, of any order, as seeds for the generation of a 
Voronoi diagram as shown in see Figure [6] 



5. Conclusions 

In this chapter statistical parameters for the truncated Pareto distribution, namely average, 
variance, median and n — th central moment have been calculated. Furthermore, also 
distribution function and survival function have been derived. These quantities allow to fit 
the various families of asteroids and the Astorb database which are characterized by a finite 
rather than infinite maximum diameter. Two possible simulations are suggested to produce 
Pareto tails. The first one results from a simple growth process, in which the increase of 
the state variable (here mass) depends on the values taken in the previous state. The second 
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one is a fragmentation process given by 3D Voronoi volumes with a fractal distribution of 
seeds. 
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